/*
**********************************************************************************
    OpenSURF642 : An TI DM642 implementation of OpenSURF.
    Copyright (C) 2011  Hai Yu

    This program is free software: you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with this program.  If not, see <http://www.gnu.org/licenses/>.
    
	This implementation is based on The original OpenSURF library created
	by Christopher Evans.The original work can be found at
		http://www.chrisevansdev.com/computer-vision-opensurf.html
	
    contact: yuhaihai12345@gmail.com	
**********************************************************************************
*/

#include "OpenSurf642.h"


void surfDetDes(    Uimage8 * image,
					IPoint ** ipts,
					Uint16  * ipts_count,
	
					Uint8     upright,
					Uint16    octaves,
					Uint16    intervals,
					Uint16    init_sample,
					ResponseDataType thres
					)
{
	IntegralImage * integral = CreateIntegralImage(image);

	FastHessians  * hessians = CreateFastHessians(	integral,
													octaves, 
													intervals, 
													init_sample, 
													thres);
	
	Surf *  surf;
	
	DetectIPoints(hessians); 
	surf = CreateSurf(integral, hessians->ipts, hessians->ipts_count);
	getDescriptors(surf, upright);
	if(ipts != 0 && ipts_count != 0)
	{
		*ipts = hessians->ipts;
		*ipts_count = hessians->ipts_count;
	}
	else
	{
		free(hessians->ipts);
	}
		
	
	DestroySurf(surf);
	DestroyFastHessians(hessians);
	DestroyIntegralImage(integral);
	
}



void surfDet(       Uimage8 * image,
					IPoint ** ipts,
					Uint16  * ipts_count,
	
					Uint16    octaves,
					Uint16    intervals,
					Uint16    init_sample,
					ResponseDataType thres
					)
{
	IntegralImage * integral = CreateIntegralImage(image);

	FastHessians  * hessians = CreateFastHessians(	integral,
													octaves, 
													intervals, 
													init_sample, 
													thres);
	
	
	
	DetectIPoints(hessians); 
	
	if(ipts != 0 && ipts_count != 0)
	{
		*ipts = hessians->ipts;
		*ipts_count = hessians->ipts_count;
	}
	else
	{
		free(hessians->ipts);
	}
		
	DestroyFastHessians(hessians);
	DestroyIntegralImage(integral);
	
}

